Cosmological constraints on non-standard inflationary quantum collapse models 
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We briefly review an important shortcoming -unearthed in previous works- of the standard version 
of the inflationary model for the emergence of the seeds of cosmic structure. We consider here 
some consequences emerging from a proposal inspired on ideas of Penrose and Diosi [JQ about a 
quantum-gravity induced reduction of the wave function, which has been put forward to address 
the shortcomings, arguing that its effect on the inflaton field is what can lead to the emergence of 
the seeds of cosmic structure [f|. 

The proposal leads to a deviation of the primordial spectrum from the scale-invariant Harrison- 
Zel'dovich one, and consequently, to a different CMB power spectrum. We perform statistical 
analyses to test two quantum collapse schemes with recent data from the CMB, including the 7-yr 
release of WMAP and the matter power spectrum measured using LRGs by the Sloan Digital Sky 
Survey. Results from the statistical analyses indicate that several collapse models are compatible 
with CMB and LRG data, and establish constraints on the free parameters of the models. The data 
put no restriction on the timescale for the collapse of the scalar field modes. 

PACS numbers: 98.80Cq,98.70Vc 



I. INTRODUCTION 

The great advances made in physical cosmology over the past few years open a window for the consideration of issues 
long dismissed as philosophical speculations. The agreement between theory and observations of the spectrum of the 
cosmic microwave background (CMB) anisotropics has strengthened the theoretical status of inflationary scenarios 
among cosmologists. The inflationary paradigm is said to account for the origin of all cosmic structure: the fluctuations 
around homogeneity and isotropy that have a fundamentally quantum nature are, according to these ideas, the seeds 
of galaxies, stars, life and humans. Hence, all structures in our universe emerge from a featureless stage described 
by a background Friedmann-Robertson- Walker (FRW) cosmology with a nearly exponential expansion driven by the 
potential of a single scalar field 1 , and from its quantum fluctuations characterized by a simple vacuum state. In 
particular, the quantum fluctuations transmute into the classical statistical fluctuations that imprint their signature 
in the CMB photons, and eventually grow into the structures we find in our universe today. 

However, when this picture is considered thoroughly, an unavoidable issue arises. According to these ideas, a 
completely homogeneous and isotropic stage, somehow evolves, after some time, into an inhomogeneous and anisotropic 
situation. Obviously, this is not simply the result of quantum unitary evolution, since, in this case, the dynamics does 
not break the initial symmetries of the system. As discussed in Ref. @ , and despite multiple claims to the contrary, 
there is no satisfactory solution to this problem within the standard physical paradigms. While much of the focus of 
the research in inflationary cosmology has been directed towards elucidation of the details of the inflationary model, 
very little attention has been given by the community to an issue of fundamental principles such as the aforementioned 
one. 

One proposal to handle this shortcoming has been developed in Refs. 0, 0. That approach attempts to deal 
with the problem by introducing a new ingredient into the inflationary account of the origin of the seeds of cosmic 
structure: the self induced collapse hypothesis, i.e. a scheme in which an internally induced collapse of the wave 
function of the inflaton field 2 is the mechanism by which inhomogeneities and anisotropies arise at each particular 



1 In the simplest models of inflation: <j>, the inflaton. 

2 As shown in [tj, one can implement the collapse hypothesis at the level of the Mukhanov Sasaki variable. 
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length scale. That work was inspired in early ideas by R. Penrose and L. Diosi [lHjj which regarded the collapse of 
the wave function as an actual physical process (instead of just an artifact of the description of Physics) and it is 
assumed to be caused by quantum aspects of gravitation. We will not recapitulate the motivations and discussion of 
the original proposal but instead refer the reader to the above mentioned papers. 

The way we treat this process is by assuming that at a certain stage in cosmic evolution there is an induced jump 
in a state describing a particular mode of the quantum field, in a manner that is similar to the quantum mechanical 
reduction of the wave function associated with a measurement, but with the difference that in our scheme no external 
measuring device or observer is called upon as triggering such collapse (as there is nothing in the situation we are 
considering that could be called upon to play such a role). The issue that then arises concerns the characteristics 
of the state into which such jump occurs. In particular, what determines the expectation values of the field and 
momentum conjugate variables for the after-collapse state. Previous works by people in our group have extensively 
discussed both the conceptual and formal aspects of that problem, and the present manuscript will not dwell further 
into those aspects, except for a very short review. 

In Ref. [5j two schemes were considered; one in which, after the collapse, both expectation values are randomly 
distributed within their respective ranges of uncertainties in the precollapsed state, and another one in which it is 
only the conjugate momentum that changes its expectation value from zero to a value in its corresponding range, as 
a result of the collapse. We will discuss the motivations and detailed characterizations of the two processes in the 
following sections. As reported in Q, the different collapse schemes give rise to different characteristic departures 
from the conventional Harrison-Zel'dovich (HZ) flat primordial spectrum. 

These aspects have been preliminarily analyzed in Ref. (Toj using a simple approach which ignores the late-time 
physics effects and simply considers what can be a reasonable expectation for the allowed deviation from the conven- 
tional flat spectrum. The main objective of this article is to consider in detail the shape of the primordial spectrum 
emerging from such a scheme, and in particular to explore with precision the deviations that should be expected once 
one takes into account the physics associated with plasma dynamics, which is responsible for the generation of the 
acoustic peaks, and other modifications to the spectrum associated with well established physics. All of these must 
be taken into account so that one can compare directly the predictions corresponding to a particular collapse model 
to recent data from the CMB fluctuation spectrum and the matter power spectrum from recent galaxy surveys, and 
thus be able to put bounds on the model's parameters. 

The paper is organized as follows. In Sect. [Til we describe the theoretical framework in which to study the models 
to be tested. Section UTTl gives some specifications about the proposed collapse schemes together with a brief physical 
motivation for them, while Sect. IIVI is dedicated to general predictions of these models, and a broad comparison to 
observational data. We define a fiducial model which is determined by the best fit values obtained by the WMAP 
collaboration with the same data used in this paper and without assuming any collapse of the wave function. It should 
be noted that the fiducial model has a value of the spectral scalar index n s different than I, while calculations we 
have carried out here of the primordial fluctuation spectrum for models including collapse are restricted to the case 
n s = 1. In Sect. [V] we use data from WMAP 7-year release and the power spectrum of the Sloan Digital Sky Survery 
DR7 LRG, to put bounds on the free parameters of the collapse models, allowing also other cosmological parameters 
to vary. In Section IVTI we discuss our results and show our conclusions. 

II. THE FORMALISM / THEORETICAL MODEL 

The starting point is the same as for the standard approaches, and in particular we will focus on one of the simplest 
inflationary models corresponding to a single scalar field, minimally coupled to gravity, with an appropriate potential. 
The action for the theory is 

S[cf>, g ab ] = J d 4 x^—g (^i^] - lv a ^V b ^ - V (4>)) . (1) 

As it is customary in such studies, we separate the fields into a "background" part, taken to be homogeneous and 
isotropic, FRW universe driven by an equally homogeneous and isotropic configuration of the inflaton field, and the 
perturbations (or "fluctuations") 3 . In this way, the metric and the scalar field are written as: g = go + Sg and 
<f> — 4>q + 8(f). This perturbative treatment requires, as usual, to deal with the gauge freedom, which we do by fixing 



In fact, the proposal has been developed to fit within a semiclassical treatment where gravitation is treated at the classical level but 
the scalar fields are completely described using quantum field theory on curved space-time in a self-consistent approach adapted to 
incorporate the collapse hypothesis. For a detailed description of this approach see 
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the gauge (conformal Newton gauge). In the present work, we will ignore the vector and tensor parts of the metric 
perturbations. The space-time metric is described as: ds 2 — a(rj) 2 [—(1 + 2<fi)dr] 2 + (1 — 2^)6^ dx i dx^\ , where *Sf and 
(j) are the Newtonian potentials. In this case, as it is well-known [B[ the field equations imply cj) = W . The scale factor 
a is normalized so that today a = 1. 

In order to further specify the setting of the problem, let us recall the standard inflationary version of cosmology. 
The universe is thought to go through the following epochs: 1) A quantum gravity regime which leads to 2) a 
preinflationary classical space-time which in turns enters 3) an inflationary regime with a = asi (BI stands for 
"beginning of inflation") and that ends when a — clei (EI stands for "end of inflation"), passing through a rapid 4) 
reheating period, taken in our approximation to occur instantaneously. The next epoch is 5) a radiation-dominated 
era starting at a = Oej, with the radiation characterized by a temperature T = Tqut, and ending at the transition 
a eq where 6) the matter (and later also the dark energy) era begins. Particularly interesting for us will be the 
recombination time occurring during the radiation-dominated regime, with a = and whose intersection with our 
past light-cone defines the surface of last scattering that is observed today in the CMB. 

We will use a different conformal time coordinate for each of the eras. Thus, we will use r\ as the conformal time 
during inflation, fj the conformal time during radiation epoch and fj the conformal time during the matter and dark 
energy epochs, and therefore the values of these coordinates at the transition points will generally differ. 

We choose the origin of rj so that during inflation we have a = jj-p^, thus inflation will end at r) = t/ei = u7a E i < ^' 
On the other hand, we choose the radiation-dominated era to start at fj — so that during the radiation epoch 
a = Cfj + B, with C = [^^-{pRadd 4 )} 1 ^ 2 and B = clei, where pR a d is the radiation energy density so that pRadO- 4 is 
constant during the era of radiation domination. Note that we are considering the reheating to occur instantaneously 
and hence the end of inflation (at rj = T)ei) and the start of radiation dominated era (at fj = 0) coincide, so that 
B = a E i- 

The study will rely on Einstein Field Equations (EFE) to zeroth and first order. The zeroth order gives rise to 
the standard solutions in the inflationary stage, where a(rj) = —jfp^i with H 2 ~ (87r/3)GV^ with the scalar field 4>o 

in slow- roll regime, so that cj)' ~ —jffj^', and the first order EFE leads to an equation relating the gravitational 
perturbation and the perturbation of the field: 

V 2 * = 4TtG<p' 8<t)' = s5(f>', (2) 

where s = 4-7tG0q. Next, we must consider the quantization of the inflaton perturbation field. It is convenient to 
work with the rescaled field y — a5<fi. In order to avoid infrared problems, we consider the system restricted to a box 
of side L, where we impose, as usual, periodic boundary conditions. We thus write the fields as 

V(V, ^>^E e 4fe %(77), 77(77, x) = -5 ^ ^'^(77), (3) 

k k 

where nk is the canonical momentum of the scaled field. The wave vectors satisfy fc^L = 2imi, with i = 1,2,3. 
Also, as usual, we write the field operators in terms of the time-independent creation and annihilation operators, 
y k (r]) = yk{v)^k + yk(v)al, and 7r fc (^) = g k {ri)a k + g k (rj)a\.. The functions yk(v), Qkiv) reflect the selection of the 
vacuum state, and here we again proceed as in the standard approaches and choose the so called Bunch-Davies 
vacuum: 

The vacuum state, defined by fife|0) = for all k, is exactly homogeneous and isotropic. We assume that at a 
certain time 77^ the part of the state characterizing the mode k (we must of course be aware that the state of the field 
is a collective state of all modes but taking here in this loose sense will do no harm) jumps to a new state, which is 
no longer homogeneous and isotropic. The detailed description of this process is as follows. 

We decompose the fields into their hermitian part as: y k = Vk(v) + wli 1 !)? ancl ^k = ^kiv) + i^Kv)- We note that 
the vacuum state |0) is characterized in part by the following: its expectation values are (y^' 1 {if)} = i^k' i 7 ))) = 0, 
and its uncertainties are Ayf J = l/2\y k \ 2 (hL 3 ) and Att^' 1 = l/2\g k \ 2 (hL 3 ). 

According to the collapse hypothesis, at some appropriate time 77S, the state undergoes an instantaneous jump to 

a different state so that the mode k is not longer in its vacuum state. 

For any state of the field |f2), we introduce the quantity d k = (f2|a^/ 7 |£l) = \d^ I \e lak so that, for that state, we 
have 

(yf) = V2M {ykd £,l ): { ^J } = ^ sr {gkd f (5 ) 
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indicating that it specifies the main quantity of interest in characterizing the state of the field. 

The analysis proceeds with the specification of the scheme of collapse determining the state of the field after the 
collapse 4 . The detailed characterization of the schemes under consideration is the main purpose of the next section. 
With such a collapse scheme at hand, one then proceeds to evaluate the perturbed metric using a semiclassical 
description of gravitation in interaction with quantum fields as reflected in the semiclassical EFE's: G a b = ^G{T a b) 
(see Eq. ([2|) for the classical first order version). To lowest order, and for the quantity of interest, this set of equations 
reduces to 

V 2 *, = s(5$' k ) a , (6) 

where (d<f)' k )Q is the expectation value of the momentum field 8<f>' k — Tt k /a[rj) on the state |0) characterizing the 
quantum part of the inflaton field. It is worthwhile emphasizing that before the collapse has occurred there are no 
metric perturbations, i.e. the r.h.s. of the last equation evaluated on the vacuum state is zero, so, it is only after 
the collapse that the gravitational perturbations appear, i.e. the collapse of each mode represents the onset of the 
inhomogeneity and anisotropy at the scale represented by the mode. Another point we must stress is that, according 
to our views, at all times, the Universe would be defined by a single state and not by an ensemble of states. The precise 
state at any time, could be written if we knew the modes that have collapsed up to that time, and the post-collapse states 
for each of these modes. The statistical aspects arise once we note that we do not measure directly and separately 
each of the modes with specific values of k, but rather the aggregate contribution of all such modes to the spherical 
harmonic decomposition of the temperature fluctuations of the celestial sphere (see below) . 

In order to be able to compare to observations we note that the quantity that is experimentally measured (for 
instance by WMAP satellite) is the anisotropy in the temperature, AT(6,(p)/T, which is expressed in terms of its 
spherical harmonic decomposition as ^2 lm ai m Yi m (9,ip). The link to theoretical calculations is made through the 
theoretical estimation of the value of the a; m 's, which are expressed in terms of the Newtonian potential on the 
2-sphere corresponding to the intersection of our past light cone with the last scattering surface (LSS): ^{t}d,Xd), 
a im = J ^( r lD,XD)Y l ^ n d 2 n. We must then consider the expression for the Newtonian potential (Eq. ©) at those 
points: 

k 

where we have introduced the factor T{k) to represent the physical effects of the period between reheating and 
decoupling. 

Writing the coordinates of the points of interest on the surface of last scattering as x = 
Rd (sin 9 sin 0, sin 9 cos (j>, cos 9), where Rd is the comoving radius to that surface, and (9, <f>) are the standard spherical 
coordinates on the sphere, and using standard results connecting Fourier and spherical expansions, we obtain 

a <™ = E TTfT I (S$' k )e in -*Y lm (9, <p)d 2 n. (8) 
k 

As indicated above, statistical considerations arise when noting that Eq. (JSJ indicates that the quantity of interest is 
in fact the result of a large number (actually infinite) of harmonic oscillators, each one contributing with a complex 
number to the sum, leading to what is in effect a two dimensional random walk whose total displacement corresponds 
to the quantity of observational interest. Note that this part of the analysis is substantially different from the 
corresponding one in the standard approach. In order to obtain a prediction, we need to find the magnitude of such 
total displacement, i.e. |a; m | 2 , and given that by assumption there are multiple random processes involved, the best 
that can be obtained is the "most likely" value for this quantity (just as in a random walk one can not expect to 
predict the actual value but at best its most likely value). We do this with the help of the imaginary ensemble of 
universes 5 and the identification of the most likely value with the ensemble mean value. We should emphasize however 
that if we knew the specific values taken by the random numbers x we will specify below, we would be able to compute 
explicitly the value of each a% m including the phase. 

As we will see, the ensemble mean value of the product (84> k )(S(j> k ')* , evaluated in the post-collapse states 6 , results 
in a form KC(k)StSn where k = TiL 3 k/{Aa 2 ) and C{k) is an adimensional function of k which encodes the relevant 



4 At this point, in fact, all we require is the specification of the expectation values of certain operators in this new quantum state. 

5 This is just a mathematical evaluation device and no assumption regarding the existence of such ensemble of universes is made or 
needed. These aspects of our discussion can be regarded as related to the so called cosmic variance problem. 

6 Note here again the difference with the standard treatment of this part of the calculation, which calls for the evaluation of the expectation 
value {S(j> k S<j>y)* on the vacuum state which as already emphasized is completely homogeneous and isotropic. 
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information of detailed aspects of the collapse scheme. Then, using the expression for |a; m | 2 , writing the sum as an 
integral, and doing a change of integration variable x = kRjj , we arrive to the following expression for the most likely 
(ML) value of the quantity of interest: 

\a lm \l 1L = A / ^Wl^fWdx. (9) 
Zira z J x 

With this expression at hand, we can compare the expectations from each of the schemes of collapse against the 
observations. We note, in considering the last equation, that the standard form of the spectrum corresponds to 
replacing the function C by a constant. In fact if one replaces C by 1 and, furthermore, one takes the function T which 
encodes the late-time physics including the plasma oscillations which are responsible for the famous acoustic peaks, and 
substitutes it by a constant, one obtains the characteristic signature of a scale invariant spectrum: |aj m |ML on m^u ■ 
Maintaining the appropriate function T leads of course to the well-known spectral shape that fits the data quite well 
once some basic cosmological parameters have been appropriately adjusted. 

In the remainder of the paper, we will focus on the effects that a nontrivial form of the function C has on the 
predicted form of the observational spectrum, and on using data to constrain aspects of the collapse models. 



III. DESCRIPTION OF THE PROPOSED COLLAPSE SCHEMES 

A collapse scheme is a recipe to characterize and select the state into which each of the modes k of the scalar field 
jumps at the corresponding time of collapse rjS. 

As we have clearly stated, we do not know exactly what kind of physical mechanism would lie behind what, at the 
semiclassical level we are working, looks like a spontaneous collapse of the wave function. Thus, specifying the collapse 
scheme should be taken as being at this point purely guesswork, which we hope to address with as much physical 
intuition as possible. The aspect that makes such efforts worthwhile is the fact that, as we shall see, the various 
specific collapse schemes lead generically to different generic patterns of deviations of the form of the primordial 
spectrum from the standard HZ scale-free shape. 

The differences in the form of the primordial spectrum lead in turn to characteristic deformations for the predictions 
of the observational spectrum which is the result of late-time and well understood physical effects on the primordial 
spectrum. Thus, observations can help us to determine which one of the naively guessed collapse schemes is favored 
by the data. 

In this paper, we focus on two schemes which we characterize simply in terms of the expectation values of the 
fundamental field operators (field and momentum conjugated) for the various modes in the state just after the 
corresponding collapse. It turns out that this limited characterization of the state (together with the times of the 
collapses) is all one needs to compute the shape of the expected spectrum. The higher level characterization of 
the states, such as the uncertainties and field-momentum correlations, turns out to be relevant only when one con- 
siders the possibility of multiple collapses per mode (see [l2j]) something that lies outside the scope of the present work. 

Scheme I: It is the scheme where both the expectation value of the field and the expectation value of the conjugated 
momentum change as a result of the collapse in a random uncorrelated manner dictated by the uncertainties of the 
precollapse state, so that, immediately after the collapse time 77^, the expectation values are determined by 

(R I) 

where x\ 2 ' are random variables, characterized by a Gaussian distribution centered at zero with dispersion equal 

to 1, Ay^ R ' IS> and Att^' 1 ^ are the uncertainties of y^' 1 ^ and n^' 1 ^ respectively in the precollapse state. This is, in 
a sense, a very simple scheme and one could argue that it represents a natural prescription that treats the field and 
conjugate momenta on an equal footing. However, a close examination (specifically, a close look at Eq. ^) indicates 
that field and momentum conjugate play rather different roles in determining the gravitational perturbation that 
results after the collapse. There we see that the momentum conjugate is the quantity that determines the Newtonian 
potential. 

This leads us to consider a collapse scheme in which it is only the expectation value of the conjugate momentum 
the one that changes as a result of the collapse (we would be taking here a view according to which, heuristically 
speaking, the uncertainty of the source of the metric perturbation, the Newtonian potential, is somehow connected 
to the triggering of the collapse of the source). See a more detailed discussion of these ideas in [5] and the motivating 
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ideas in [1 4]. Thus we define: 

Scheme II: It is the scheme where it is only the expectation value of the conjugated momentum the one that 
changes as a result of the collapse, and it does so in a random way, dictated by the corresponding uncertainties of the 
precollapse state, so that, immediately after the collapse time rj^, the expectation values are determined by 

(vl R ' I] fo£)) n = 0, (*f' J) (tf)) n = ^VK^< ™ 

( R I) 

where 

2 £ixe 

random variables, characterized by a Gaussian distribution centered at zero with dispersion of 1, and 

AttI^ is the uncertainty of fl^' in the precollapse state. 

We should note that an additional and crucial piece of information would be required to determine the primordial 
spectra: the exact values of the collapse times for each mode. Again, having no knowledge of the physics behind the 
collapse, all we can do is to assume that whatever that is, it has no intrinsic preferential directionality and thus that 
77S is only a function of k = \\k\\, but otherwise take that function as unknown. We will then parametrize in simple 
ways our ignorance about this function and attempt to extract constraints from the data. 

There are of course many more possibilities of collapse schemes but we limit our consideration to only those two, 
because they seem simple and quite natural. Next, we investigate these two schemes in detail at a quantitative level. 



IV. COMPARING WITH OBSERVATIONAL DATA 

Given the schemes proposed, and after some lengthy algebraic manipulations, one can determine (see ref. [13]) the 
functions C(fc)'s and then use them to compare the resulting predictions for the shape of the spectrum to observations. 
It was shown in Refs. [H,[l(| that we can recover the power spectrum of initial fluctuations of the standard cosmological 
model if C(k) = 1. Therefore, the power spectrum of the collapse models can be written as: 



P{k) = A s kC{k) (12) 
where A s is the amplitude of the scalar fluctuations. 

The function C(k) resulting from Eq. (flO|) for scheme I has the following form: 

Ci{k) = l + A sm 2 Afe + — sin(2A fc ), (13) 

z k z k 

whereas the one corresponding to the scheme II is: 

C u (k) = 1 + (l - \) sin 2 A fe - — sin(2A fc ), (14) 
V z k / Zk 

where Zk = krj^., Ak = k(r]Ei — ilk) an d Vei refers to the conformal time at the end of the inflationary period 7 . We 
refer to them hereafter as Model I and Model II, respectively. 

In Ref. [lOj a preliminary study of these two schemes and one additional scheme was performed. That was a 
relatively simple analysis concentrating on the main features of the resulting spectrum, but ignoring the late-time 
physics corresponding to the effects of reheating and acoustic oscillations (represented by T(k)). Thus, the actual 
comparison with empirical data was not possible except for obtaining order-of-magnitude estimates. A detailed 
comparison with the very precise data available, requires a much more complex analysis, such as the one we will 
conduct in this work. 

We recall that the standard form of the predicted spectrum is recovered if one replaces C(k) by a constant. We 
will explore the sensitivity for small deviations of the "zk independent of k pattern" by considering a linear departure 
from the situation in which Zk is independent of k. That will be characterized by Zk as Zk = A + Bk. This will allow 
us to examine the robustness of the collapse scheme in as far as predicting the standard spectrum. In such way, the 



7 This expression is different from the one in previous papers [5I llOj. because in those, the work relied in an aproximation in which the 
effects of the plasma physics relevant to the era between the end of inflation and decoupling are ignored. 
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times of collapse for each mode k can also be written in terms of A and B as follows: r)Z = A/k + B. We can use this 
formula to compute the collapse time for the relevant modes we observe in the CMB, namely those that cover the 
range of the multipoles of interest, 1 < I < 2600. We can made use of the approximate relation 8 I ~ irkRu, where 
Rd is the comoving radius of the last scattering surface, in order to interpret heuristically the result of the analysis 
and to set reasonable viability constraints on the parameters of the model. In standard cosmology 9 it is given by 
Rd = -J^ (l — y/aoj , where we have normalized the scale factor to be ao = 1 today, so ar> = a(rjn) — 10~ 3 , and Hq 
is the Hubble constant today. Its numerical value is Rd = 5816.31 hr 1 Mpc. Thus, the relevant modes for the CMB 
are those in the range 4 x 10~ 5 Mpc -1 < k < 0.11 Mpc -1 . The collapse times for these modes can be regarded as 
the times in which inhomogeneities and anisotropies first emerged at the corresponding scales. 



A. Physical meaning of the obtained results 

The analysis carried out in previous works on this approach jfj [lol ] does not take into account the evolution of the 
perturbation beyond the end of inflation and thus their results, even though they provide certain general qualitative 
and quantitative information about the modifications the general approach leads to, cannot be considered as leading to 
actual bounds on the parameters extracted from data. In this work we want to extract actual bounds that could, in the 
future, be taken as clues on the nature of the physics that might lie behind what we represent at the phenomenological 
level by the collapse mechanism. In particular, we need to ensure that when considering the collapse as being described 
within the inflationary regime, the times where the relevant collapses occur do indeed fall within that regime. Hence, 
we will impose a priori limits on the values of A and B, so that the collapse of the wave function can occur only at 
a time between the beginning (tsi) and the end (£ej) of the inflationary epoch. In terms of conformal times: 

r]Bi < V < r/Ei- (15) 

The period of inflationary expansion prior to the radiation dominated era, corresponds to negative conformal time. 
The initial singularity (or more precisely, the quantum gravity regime) is pushed back into large and negative values 
of the conformal time and can be pushed arbitrary far depending on the duration of inflation [13J . 

Now, we estimate the time scales for the beginning and the end of the inflationary epoch. The value of the 
conformal time at the end of inflation, is determined from the " temperature" of the radiation era viewed as a function 
of the scale factor: a(to)T(to) = a(tEi)T(tEi) where to is the present time, a(to) = 1, t B i is the time at the end of 
inflation, and To = 2.728A". Taking the usual values for inflationary models, T(isj) — 10 15 GeV, and consequently 
d(tEi) = 2.35 x 10 - 28 . We then use the relation between conformal time and scale factor during the inflationary era, 
a ( f l) — — TTTtj ' wnere Hi is the value of the Hubble constant during this epoch. In most inflationary models, Hj can 
be estimated as: Hj = ^pVj where V/ is the inflationary potential. Let us assume that Vj ~ 10 60 GeV 4 , then we 
obtain Hi = 2.37 x 10 11 GcV and thus 

T)ei = -1.8 x 10 16 GeV -1 = -7.2 x 10~ 22 Mpc. (16) 

We express the results in Mpc because these are the common units for time and length used in Boltzmann codes 
that solves the recombination equations, such as CAMB [lij]. In order to solve the "horizon" and "flatness" problem, 
inflationary models require 80 e-folds of inflation: 

log ( = Hifa* - ^ = 80 - 
\a{t B i)J 

Thus, we obtain: a(t B i) = 4.24 x 10~ 63 , and using the relation between the scale factor and conformal time during 
inflation, we get: 

t]bi = -9.9 x 10 50 GeV -1 = -1.6 x 10 14 Mpc. (18) 



The relation between the angular scale 8 and the multipole I is 8 ~ it/I. The comoving angular distance, c2a, from us to an object of 
physical linear size L, is d\ = L/(a8). L/a ~ 1/fc, djy = Hp if the object is in the LSS, and using the first expression in this footnote, 
we get I ~ irkRjj. 

In the present work we will be ignoring the effects of the late time acceleration associated with the so called "Dark Energy" as this 
complication is thought not to impact on the results in any substantial way. 
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Let us recall now that in the models that are being studied, each mode is thought to collapse at a different time 
given by 

Vk = j + B, (19) 

where A and B are constants. Therefore, for the collapse of the wave function to occur during the inflationary period, 
the values of A and B should satisfy the following relation for all values of k: 

A A 
Vbi - j < B <r] EI - -. (20) 

Demanding that these inequalities hold for all the relevant modes, leads to the desired a priori bounds on the model 
parameters. 



B. Characteristic signatures of the models on the CMB fluctuation spectrum 



In order to analyze the effects on the CMB fluctuations power spectrum, let us first define the fiducial model, which 
will be taken just as a reference to discuss the results we obtain for the collapse models. The fiducial model is a ACDM 
model with the following cosmological parameters: baryon density in units of the critical density, Vlsh 2 = 0.02247; 
dark matter density in units of the critical density, QcDAih 2 = 0.1161; Hubble constant in units of Mpc _1 km s _1 , 
Hq = 68.7; reionization optical depth, t = 0.088; and the scalar spectral index, 0.959. These are the best-fit 

values presented by the WMAP collaboration using the final 7-year release data [15[ and the power spectrum from 
Sloan Digital Sky Survery DR7 LRG [3. 




FIG. 1: The temperature autocorrelation (TT) power spectrum for Model I (Left: A — —10, Right: A — — 10 8 ). All models 
are normalized to the maximum of the first peak of the fiducial model. The value of \ 2 is calculated using only WMAP 7- year 
release data. The solid line corresponds to the fiducial model. 

Figures Q] [5] and [6] show the prediction for the temperature fluctuation spectrum for Models I and II for various 
values of the parameters A and B. The value of % 2 in the figures is calculated using both temperature and temperature- 
polarization fluctuation data from WMAP 7-year release. We can see that the different collapse models have distinct 
effects on the CMB temperature fluctuation spectrum, and furthermore, that these effects depend rather strongly on 
the values of the parameters A and B which determine the collapse times of the various modes. 

Let us first analyze the case of Model I: When fixing A — —10 or A = — 10~ 4 , we find that there is no change within 
the range of B studied here (with respect to the fiducial model) in the position of the Doppler peaks or in the height 
of the first peak. On the other hand, we note an increase in the height of the secondary peaks, with the amplitude 
depending on the value of the free parameter B. Next, we set A = — 10 s (see Fig.[T] right panel); in this case there is 
no change in the first peak (with respect to the fiducial model), while there is a extremely small increase in the height 
of the second peak and a a similar decrease in the height of the third peak. In fact one does not observe any important 
change in the spectrum for different values of B (other values of B in the explored range give the same spectrum). 
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FIG. 2: The temperature autocorrelation (TT) power spectrum for Model I (A — — 10~ 4 ). All models are normalized to the 
maximum of the first peak of the fiducial model. The value of \ 2 is calculated using only WMAP 7-year release data (both 
temperature and temperature-polarization power spectrum are included). The solid line corresponds to the fiducial model. 
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FIG. 3: The temperature-polarization (TE) cross-power spectrum for Model I (Left: A — —10, Right: A = — 10 8 ). All models 
are normalized to the maximum of the first peak of the fiducial model. The value of \ 2 is calculated using only WMAP 7- year 
release data (both temperature and temperature-polarization power spectrum are included). The solid line corresponds to the 
fiducial model. 



The small difference between the spectrum of the collapse models and the fiducial model is due to the fact that the 
former uses the canonical value for the spectral index (n s = 1) while the latter is the best fit value to the data found 
by the WMAP collaboration assuming the standard cosmological model (which differs slightly from 1). We have also 
studied the results for other various large negative values of A (not shown here) and found that the behavior is similar 
to that of the case A = — 10 8 . This is the result of a simple preliminary analysis; the exact value of A at which 
the behavior of the CMB spectrum changes becomes important, phenomenologically speaking, will be determined 
from the statistical analysis (see Section [Vj. As the most accurate observational data corresponds to the first peak, 
it is expected that the collapse Model I with large values of A will be phenomenologically indistinguishable from the 
fiducial model, as seen in Fig. Q] (right panel). Indeed, the difference in the value of the quality-of-fit estimator x 2 , is 
found to be not significant within la. 

Figs. [3] and UH show the prediction for the temperature-polarization cross-power spectrum for Model I. Here again, 
the models with A = —10 s differ very little from the fiducial model and show no relevant dependence of the results on 
values of B within the explored range. Regarding the cases of A = —10 and A = — 10 -4 , we find (in comparison with 




FIG. 4: The temperature-polarization (TE) cross-power spectrum for Model I (A — — 1CP 4 ). All models are normalized to the 
maximum of the first peak of the fiducial model. The value of \ 2 is calculated using only WMAP 7-year release data (both 
temperature and temperature-polarization power spectrum are included). The solid line corresponds to the fiducial model. 



the fiducial model) an increase in the values of the secondary peaks and a decrease in the amplitude of the valleys, 
but the relative change is less relevant than that in the temperature fluctuation power spectrum. 

Let us now focus on Model II (see Fig. [SJ. Here again, the effect is different for different values of A: Considering 
firstly the case of A = —10, we find that there is a shift in the position of the peaks, and in one case (B = — 10 3 ), the 
first peak is replaced by two peaks. On the other hand, when setting A — — 10 8 or A = — 10~ 4 , we find generically 
a shift in the position of the first peak. In both cases, there is also a change in the height of the secondary peaks, 
with the magnitude of the change depending on the values of A and B. Figs. [7] and [8] show the prediction for the 
temperature-polarization cross-power spectrum for Model II. Here again, we observe (in comparison with the fiducial 
model) an increase in the value of the peaks and a decrease in the values at the valleys for all cases, with the magnitude 
of the changes depending on the value of B. Therefore, in this case, collapse models can be clearly distinguished from 
the fiducial model. The difference in behavior of the two models can be understood by looking at Eqs. (fT3| and (|14jl : 
for sufficiently large values of A and/or B, Ci(k) becomes independent from k, while the second term of C%(k) does 
not vanish in the limit of large values of hrj^ . We would like to emphasize the importance of this preliminary analysis 
on the behavior of the collapse schemes, to determine the appropriate method for the statistical analysis to explore 
the various relevant regions of parameter space, as we described in detail in Sect. [Vj 



V. RESULTS 



The observational data used for the analysis are the temperature and temperature-polarization power spectra 
obtained from the final WMAP 7-year release fl5ll. and other CMB experiments such as CBI [17| , ACB AR pi" 



BOOMERANG [I9j,l2fj], BICEP [2l| and QUAD 22], together with the matter power spectrum traced by LRGs as 
measured in the Sloan Digital Sky Survey DR7 [lf|. We consider a spatially-flat cosmological model with adiabatic 
density fluctuations, in which we add the effect of the collapse models in the power spectrum of the initial fluctuations. 
The parameters allowed to vary are: 

p = (n B h 2 , n CD Mh 2 , e, T , a„a, b) , (21) 

where O is the ratio of the comoving sound horizon at decoupling to the angular diameter distance to the surface of 
last scattering, r is the reionization optical depth, A s is the amplitude of the primordial density fluctuations, and A 
and B are the model parameters related to the conformal time of collapse of each mode (see Eq. (|T9|) ). Given that the 
primordial power spectrum of density fluctuations for the collapse models was computed only for the scale invariant 
case [1, [nj, the scalar spectral index of density fluctuations (n s ) is fixed to 1 in this paper. However, we should keep 
in mind that the collapse models allow different values for n s and the corresponding power spectrum will be studied 
in future works. 

In order to place constraints on the parameters of the quantum collapse models, we modified the primordial power 
spectrum according to the schemes described in Section [II] We performed our statistical analysis by exploring the 




FIG. 5: The temperature autocorrelation (TT) power spectrum for Model II (Left: A = —10, Right: A = —10 s ). All models 
are normalized to the maximum of the first peak of the fiducial model. The value of \ 2 is calculated using only WMAP 7-year 
release data (both temperature and temperature-polarization power spectrum are included). The solid line corresponds to the 
fiducial model. 
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FIG. 6: The temperature autocorrelation (TT) power spectrum for Model II (A = — 1CP 4 ). All models are normalized to the 
maximum of the first peak of the fiducial model. The value of x 2 is calculated using only WMAP 7-year release data (both 
temperature and temperature-polarization power spectrum are included). The solid line corresponds to the fiducial model. 



parameter space with Monte Carlo Markov chains generated with the publicly available CosmoMC code of Ref. [23[ 
which uses the Boltzmann code CAMB [l4[ to compute the CMB power spectra. 

In a first trial, we performed the statistical analysis without imposing any prior on the baryon density and found 
that the confidence interval obtained for fi^/i 2 turned out to be very large compared to the one obtained by the 
WMAP collaboration (i.e. without considering the collapse scheme). Therefore, we introduced a gaussian prior on 
the baryon density, using an independent data, namely the Big Bang Nucleosynthesis (BBN) bound pil . [25| . 

Initially, we intended to perform the statistical analysis allowing the two parameters of the model, A and B, to 
vary independently. However, after the first runs we realized that for a very large range of values of A, there was 
always a range of values of B providing a good fit to the data. Furthermore, in order to allow the Markov chains to 
explore various orders of magnitude in those parameters, we used the following reparametrization: a = sinh(A) and 
b = sinh(_B) (another possible reparametrization is a — \og(A) and b = log(_B), but we are interested in exploring both 
positive and negative values of A and B). We found that the Markov chains did not converge when both a and b were 
allowed to vary independently. The reason for this is simply the fact that there are several maxima of the probability 
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FIG. 7: The temperature-polarization (TE) cross-power spectrum for Model II (Left: A — —10, Right: A — — 10 8 ). All models 
are normalized to the maximum of the first peak of the fiducial model. The value of \ 2 is calculated using only WMAP 7-year 
release data (both temperature and temperature-polarization power spectrum are included). The solid line corresponds to the 
fiducial model. 
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FIG. 8: The temperature-polarization (TE) cross-power spectrum for Model II (A — — 10 -4 ). All models are normalized to the 
maximum of the first peak of the fiducial model. The value of x 2 is calculated using only WMAP 7-year release data (both 
temperature and temperature-polarization power spectrum are included). 



function in the parameter space, but the method of Markov chains does not allow to explore all maxima at the same 
time. In contrast, we found that when we fixed values of A within a rather large range of values, convergence of the 
Markov chains is found. Moreover, let us recall that the primordial power spectrum in collapse models is a function 
of Zk = A + Bk; by taking B = 0, we recover the spectrum corresponding to the standard ACDM model with spectral 
index n s — 1, up to an overall normalization factor. Thus, the parameter B gives a measure of the departure of the 
primordial power spectrum with respect to that of the standard model, while phenomenologically A behaves as a 
normalization factor that is degenerated with A s . We should note, however, that although the effect on the spectrum 
of A and A s is very closely connected, their physical interpretation is quite different. The results of the statistical 
analysis, for fixed values of A are shown in Tables HI [TT1 and IPLTl for collapse models I and II. 

Within Model I, a preliminary analysis of the modifications of the CMB fluctuation power spectrum associated 
with the collapse models indicates that for large values of the parameter A, there is no dependence on the value of B, 
and also that there is a very small difference between those models and the fiducial one. This last feature is related 
to the fact that the fiducial model involves a different value of the spectral index n s (see discussion in Section IIVBI) . 
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TABLE I: Model I. Mean values and la error for b = sinh B considering different fixed values of parameter A. The x 2 estimator 
is computed for the whole data set. For the fiducial model, Xmin = 3909. Column 4 refers to the viability of the model (collapse 
times happen during inflation), and column 5 shows whether there is any restriction to the viability in the range shown by 
column 2. 
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These results are, indeed, reflected in the statistical analyses by the fact that for |^4| < 20, nontrivial bounds on B 
can be found, while for larger values of A, any value of B gives a good fit to the data. Furthermore, the values of 
the corresponding cosmological parameters obtained from the statistical analyses for models with \A\ > 20 is within 
the la bounds established by the standard analysis of WMAP collaboration made without considering the collapse 
scheme. In order to confirm the persistence of this behavior, we have explored a large range of A values, from — 10 9 
to 10 9 , changing this value in each step by one order of magnitude, and found results consistent with the previous 
conclusion. In Table U we show the bounds on B obtained for fixed values of \A\ < 20 for Model I. 

A similar analysis performed for Model II, shows that we should not expect the same behavior as the one obtained 
for Model I. Indeed, Fig. [5] shows (for two different values of A) that the prediction of Model II on the GYs depends 
rather strongly on the value of B. Results from the statistical analysis for fixed values of A are shown in Table HT1 and 

nn 

As mentioned previously in this section, we could not perform a statistical analysis with A and B varying jointly as 
free parameters, due to the fact that there are several maxima of the probability function. Indeed, Tables [TT1 and IIIII 
show that, for several values of A, we obtain two different maxima for which the difference between the value of Xmin 
does not exceed the value of 1 and thus one cannot discriminate between them. There is even one case (A = 0.1) 
in Model II (see Table ITTTT) in which there are three maxima. However, we should emphasize an important difference 
between the statistical analysis performed for A = —300 in Model II and the ones corresponding to the other values 
mentioned above. In the first case, we have performed two different statistical analyses with different initial values 
of b = sinh(_B) when running COSMOMC and obtained convergence of the Markov chains to one maximum at each 
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TABLE II: Model II. Mean values and la error for b — sinhB considering different fixed negative values of parameter A. \ 2 is 
calculated for the whole data set. For the fiducial model, Xmin = 3909. Column 4 refers to the viability of the model (collapse 
times happen during inflation), and column 5 shows whether there is any restriction to the viability in the range shown by 
column 2. 
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time. The values of the errors have been calculated with the GETDIST program, as is usually done in this kind of 
analyses. On the other hand, for values of \A\ < 1, we have performed one statistical analysis for each value of A and 
obtained a marginalized likelihood with two peaks, for the parameter B. We have tried to change the initial value of 
b = sinh(B) in order to obtain convergence to a single value at the time, with no success. As the GETDIST program 
determines the confidence interval assuming that there is just one maximum of the likelihood function, we need some 
alternative method to calculate those. We have thus computed the confidence interval for values of \A\ < 1 in Model 
II by estimating the limits of integration for which the integral of the marginalized likelihood function yields 68% of 
the value of the same integral over that peak. Therefore, it is not surprising that the la errors calculated for A < 1 
are 4 times larger than those calculated by the GETDIST program. 

Next, let us recall the discussion of Sect. IIV Al where for a given value of A, we calculated the corresponding range 
of values of B that ensure that the time of collapse for the mode k occurs during the inflationary period. Then, by 
considering the values of A for which the collapse schemes are tested (see Tables HI [TT] and llllj) , Eq. (PO)) implies that 
for A > we must have B < and that for A < the condition is B < —j^— .Therefore, we have added the line 

corresponding to B = in Figs. [9] and [TTIand the line corresponding to B = in Fig. [T0l and [T2l above 
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TABLE III: Model II. Mean values and la error for b — sinhB considering different fixed positive values of parameter A. \ 2 is 
calculated for the whole data set. For the fiducial model, Xmin = 3909. Column 4 refers to the viability of the model (collapse 
times happen during inflation), and column 5 shows whether there is any restriction to the viability in the range shown by 
column 2. 
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which the solutions found are nonphysical and thus excluded. It should be noted that there are some cases where only 
parts of the solution are physically viable. In order to facilitate this discussion we have added column 4 in Tables HI 
III1 and Hill which indicates if the model is viable or not, and column 5 in the same tables indicates the cases where 
there is a further restriction on the range of values of B for the viability of the solution. In the rest of this section 
we will discuss only those results that are relevant for our model, and exclude those nonviable values mentioned in 
the above paragraph. Recall that, unfortunately, and as it has been already mentioned, the present analysis does 
not allow us to determine any possible degeneration, as far as the fit is concerned, between the model parameters A 
and B. However, from Fig. [9] we can see that the allowed values of B - within la error- for A > —20 turn out to lie 
in the range | sinh(B)| < 3.2 for Model I; recall that for A < —20, the resulting CMB spectrum is the same for all 
values of B and very similar to the fiducial model and therefore, any value of B among those tested in this paper 
(B = — 10 9 • • • 10 9 )) provides a good fit to the data. On the other hand, for Model II, the allowed values of B - within 
Icr error- lie in the range | sinh(_B)| < 3 for Model II for all values of A studied in this paper. This allows us to set 
the bounds B < 1.88 Mpc for A < -20 in Model I, and B < 1.81 Mpc for Model II. On the other hand, it is also 
interesting to note that for all the studied cases with \A\ < 1 in Model I, the confidence interval includes the case 
B = with a maximal allowed departure from that value of 0.6. In contrast, for the cases \A\ > 1 in Model I and 
for almost all cases considered in Model II, the value B = is excluded within 2a, except for the case A = — 10 5 for 
Model II and A — 10 for Model I. 

The results of the statistical analysis for the cosmological parameters are shown in Figs. [T3l[T4l [T5l andlTBlfor Models 
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I and II. We can distinguish two behaviors: i) Model I with —20 < A < — 1; Model II with A < — 1; ii) Models I with 
\A\ < 1; Model II with \A\ < 1. Let us remind that for |A| > 20 in Model I we recover the scale-invariant HZ spectrum 
and therefore the constraints on the cosmological parameters are those estimated by the WMAP collaboration (lij . 
For models included in case i), most of the estimated values for the cosmological parameters within the collapse models 
are in agreement with those obtained by the WMAP collaboration using the 7-year data release and considering a 
standard ACDM model. The values obtained for 8 are marginally consistent with those obtained by the WMAP 
collaboration. However, there is a better agreement in the results obtained for Hq, which is derived from and 
other cosmological parameters. Let us now discuss some exceptions. For the case A = —10 in Model I, the values 
obtained for Vt^h 2 and Hq are at variance with the values obtained by the WMAP collaboration within 2a, while 
there is agreement within 3a. For the case A = — 10 5 in Model II, the values obtained for fi^/i 2 and H are at 
variance with the values obtained by the WMAP collaboration within la, while there is an agreement within 2a. 
On the other hand, it is important to discuss consistency of our results with independent data such as the baryon 
density inferred from BBN [24], 10 and the constraint on the Hubble constant presented in Ref. [26|. It should be 
noted that the la BBN region for f^/i 2 shown in Figs. H"3 l IT4" 1 [T5l and [TBI refers to the constraint obtained using the 
observational abundances of deuterium, which is the most stringent one. However, one should keep in mind that by 
considering the 4 He abundance, the constraint on the baryon density softens substantially, and in fact it leads to a 
region which is consistent with the constraint obtained with the CMB data. Indeed, the values we obtained for il^h 2 
are marginally consistent with the value inferred from deuterium abundance BBN constraints within la, while there 
is full consistency with value inferred from the primordial abundance of He. On the other hand, the values of Hq 
obtained in our study are marginally consistent with the constraint extracted from Hubble Space Telescope (HST) 
data. 

Let us now discuss models included in case ii). The values of Qi>h 2 obtained are consistent with the values obtained 
by the WMAP collaboration in the context of the standard cosmological model within 2a. On the other hand, there 
is disagreement with the deuterium inferred BBN constraint on Qi,h 2 within 3a, while there is consistency with the 
value inferred from the 4 He abundance. The values of Hq obtained are at variance with the ones obtained by the 
WMAP collaboration within la, while there is agreement within 2a and with the constraints obtained with the HST 
data. The values of O also differ at la. Other values of cosmological parameters obtained do agree with the ones 
obtained by the WMAP collaboration in the context of the standard inflationary scenario. On the other hand, for 
both models, we found that there is no agreement with the values of A s obtained in almost all of the cases of fixed 
values of A explored. This can be understood as due to the fact that the effective amplitude of the power spectrum 
of density fluctuations depends on A and A s , so the parameters act at a certain level as degenerated parameters. 
However, there is agreement between the values for erg 11 obtained and those of the fiducial model. Note, however, 
that this is a derived parameter which depends on the effective amplitude. It should also be noted that the models 
included in case i) are preferred by the data as compared to the models included in case ii), as can be inferred from 
the value of Xmin °f the respective statistical analyses listed in Tables HI |H] and IIIII 
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FIG. 9: Results for Model I: Bounds on b — sinh(B) obtained for fixed negative values of A. The area below the (blue) dotted 
line sinh(i?) < sinh(— k A ) indicates the region where the collapse of all modes happens during the inflationary period. 



10 Recall that have put a gaussian prior on the baryon density using constraints inferred from BBN. 

11 as is the rms mass fluctuation amplitude in spheres of size 8 h _1 Mpc. 
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FIG. 10: Results for Model I: Bounds on b = sinh(i?) obtained for fixed positive values of A. The area below the (blue) dotted 
line sinh(_B) < sinh(— k A ) indicates the region where the collapse of all modes happens during the inflationary period. 
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FIG. 11: Results for Model II: Bounds on b — sinh(B) obtained for fixed negative values of A. The area below the (blue) dotted 
line sinh(B) < sinh( ~ A ) indicates the region where the collapse of all modes happens during the inflationary period. 

When analyzing the correlations between the cosmological parameters flbh 2 , Q c h 2 and 8, r and B, we find that for 
both Models I and II the values of the correlation coefficient pse, PBn b h 2 an d Pbh are positive for positive best- fit 
values of B, while pssi^h 2 is positive for negative best-fit values of B. Furthermore, the correlation is strong for Vt^h 2 , 
il c h 2 and Hq while for 9 and r the correlation coefficient is less than 0.3 for all the cases studied in this paper. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the phenomenological predictions of the collapse models developed in Refs. 0,113, H3] 
and incorporated the effect of late-time plasma physics. We have calculated the prediction for the CMB fluctuation 
spectrum in the case where the collapse of the k modes of the scalar field during the inflationary period is included. 
Furthermore, we have performed a statistical analysis in order to compare the prediction of the collapse models with 
recent data from the CMB fluctuation spectrum and the matter power spectrum obtained by the SDSS collaboration, 
setting bounds on the model parameters A and B that characterize the collapse times of the scalar field modes. 
Results from the statistical analyses were discussed in Sect. IVland can be summarized as follows: 

• For \A\ > 20 in Model I, any value of B provides a good fit to the data within the range studied in this paper 
(B = -10 9 ---10 9 ). 
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FIG. 12: Results for Model II: Bounds on b = sinhB obtained for fixed positive values of A. The area below the (blue) dotted 
line sinh(B) < sinh(— k A ) indicates the region where the collapse of all modes happens during the inflationary period. 



• For \ A\ < 20 in Model I, there is a nontrivial range of values of B that provide a good fit to the data. However, 
values of A lying in the range: 10~ 4 < A < 20 are not viable since the collapse of the modes occurs after the 
inflationary period. Furthermore, the data indicate that the values in the range —20 < A < — 1 are preferred. 

• The best fit values of B obtained for |A| < 20 in Model I lie in the range B < 1.88 Mpc. 

• In Model II, there is a nontrivial range of values of B that provide a good fit to the data for all values of A 
tested (A = — 10 9 • • • 10 9 ). However, values of A lying in the range: 10~ 4 < A < 10 9 are not viable since the 
collapse of the modes occurs after the inflationary period. Furthermore, the data indicate that the values in the 
range \A\ > 1 are preferred. 

• The best fit values of B obtained for all values of A studied in this paper for Model II lie in the range B < 
1.88 Mpc. 

• The value obtained for the cosmological parameters is consistent within 3cr with those obtained by the WMAP 
collaboration assuming a standard inflationary scenario and also with bounds established by BBN and the value 
of H obtained with the HST. 

This analysis allows us to compare the value of the scale factor at the collapse time cl(j]^,), with the traditional 
value of the scale factor at "horizon crossing" which is often set to mark the quantum to classical transition in the 
standard explanation of inflation: a?. The "horizon crossing" occurs when the length corresponding to the mode k 
has the same size that the Hubble radius HJ 1 (in comoving modes, k = aHi) therefore, ajf = a(i]J?) = -gU = g ^ v . 
Thus the ratio of the value of the scale factor at horizon crossing for mode k and its value at collapse time is 

= krgVk) = A + Bk. (22) 

a l 

The results discussed in Sect. |V]show that two cases are consistent with the data: the one for which the collapse 
time of inflaton field modes is previous to the time of "horizon crossing" of all modes of the inflaton field, and the 
opposite case. However, the data seem to favor the case in which the collapse happens before the modes cross the 
horizon. 

As we found that there are good fits (for suitable ranges of B) for all values of the parameter A, we must conclude 
that the existing data do not lead to a preferred scale for the times of collapse of the wave function for the relevant 
modes. However, for each value of A, interesting bounds can be placed on the value of the parameter B which in this 
approach characterizes the modifications of the primordial spectral shape with respect to the HZ conventional flat 
scale-free spectrum. 
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FIG. 13: Results for Model I: Best-fit parameter values and la errors for the cosmological parameters obtained for fixed negative 
values of A. Comparison with the values obtained by the WMAP collaboration (and other data sets, where relevant), is also 
shown. For convenience, the value of —A is plotted in the x-axis. Values of —20 < A < — 1 belong to case i) described in the 
text, while values of \A\ < 1 correspond to case ii). 



Acknowledgments 

The authors are grateful to Gabriel Leon for useful discussions about the collapse models. Numerical calculations 
were performed with the KANBALAM facility located at UNAM. The authors would like to thank the people of 
DGSCA-UNAM for computational and technical support. S. J. L. is supported by PICT 2007-02184 from Agenda 



20 



0.025 



G 0.022 



76 



CMB + LRG 1-0 region 
BBN 1-c region 



l [ 



1(T 
Value of A 



CMB + LRG 1-c region 
HST 1-c region 



10" 
Value of A 



0.124 
0.122 

0.12 
0.118 
0.116 
0.114 
0.112 

0.11 
0.108 





CMB + LRG 1 -c region 























10"° 
Value of A 



CMB + LRG 1-o region 



10"' 
Value of A 



0.12 



0.1 



CMB + LRG 1-c region 



10" 
Value of A 



0.88 r 
0.87 - 
0.86 - 
0.85 - 
0.84 - 
0.83 - 
0.82 - 
0.81 - 
0.8 - 
0.79 L*. 



CMB + LRG 1-o region 



10" 
Value of A 



FIG. 14: Results for Model I: Best-fit parameter values and la errors for the cosmological parameters obtained for fixed positive 
values of A. Comparison with the values obtained by the WMAP collaboration (and other data sets, where relevant), is also 
shown. All of these values correspond to case ii) described in the text. 
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FIG. 15: Results for Model II: Best-fit parameter values and 1<t errors for the cosmological parameters obtained for fixed 
negative values of A. Comparison with the values obtained by the WMAP collaboration (and other data sets, where relevant), 
is also shown. For convenience, the value of —A is plotted in the a>axis. Values of \A\ > 1 belong to case i) described in the 
text, while values of \A\ < 1 correspond to case ii). 
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FIG. 16: Results for Model II: Best-fit parameter values and 1<t errors for the cosmological parameters obtained for fixed 
positive values of A. Comparison with the values obtained by the WMAP collaboration (and other data sets, where relevant), 
is also shown. All of these values correspond to case ii) described in the text. 
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